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Abstract. This paper presents results from modelling 450 p.m and 850 fim continuum and HCO + line observations of three 
distinct cores of an infrared dark cloud (IRDC) directed toward the W51 GMC. In the sub-mm continuum these cores appear 
as bright, isolated emission features. One of them coincides with the peak of 8.3 yum extinction as measured by the Midcourse 
Space Experiment satellite. Detailed radiative transfer codes are applied to constrain the cores' physical conditions to address 
the key question: Do these IRDC-cores harbour luminous sources? The results of the continuum model, expressed in the \ 2 
quality-of-fit parameter, are also constrained by the absence of 100 [im emission from IRAS. For the sub-mm emission peaks 
this shows that sources of ~ 300 Lq are embedded within the cores. For the extinction peak, the combination of continuum and 
HCO + line modelling indicates that a heating source is present as well. Furthermore, the line model provides constraints on 
the dumpiness of the medium. All three cores have similar masses of about 70 - 150 Mq and similar density structures. The 
extinction peak differs from the other two cores by hosting a much weaker heating source, and the sub-mm emission core at the 
edge of the IRDC deviates from the other cores by a higher internal dumpiness. 

Key words. ISM: clouds - ISM: dust, extinction - ISM: structure - Stars: formation - Submillimeter 



1. Introduction 

InfraRed Dark Clouds (IRDCs) have b een observed in re- 
cent mid-IR surve ys of the Galactic Plane JPerault et all 1 19961 
lEean et all Il998h as prominent absorption features against 
the bright mid-IR Galactic emission. Due to lack of 60 fim 
and 100 urn emis sion, IRDCs must be cold, i.e., below 13 K 
fegan et allfl9 98). This cold nature was confirmed, by follow- 
up observations of molecular lines and sub-mm continuu m ob- 
serva t ions and placed mos t IRDCs below 15 K dCarev e t all 
Il998t iTevssier et all 12002). Furthermore, these line observa- 
tions provided radial velocities and hence estimates of the dis- 
tances. Most IRDCs are at distances above 1 kpc and are pref- 
erentially l ocated toward spiral arms and the Galactic Ring 
fegan et allll999tlCarev et alll2000al) . 

Thousands of IRDCs have been found throughout the 
Galactic Plan e and counterp a rts are found at diffe rent wave- 
lengths JSimon et all 1200 lb . ICarev et all feOOObl) identified 
a number of IRDCs from th e Midcourse Space experiment 
(MSX) Galactic Plane survey JPrice et alll200ll) and observed 
these IRDCs in 850 yum dust continuum emission. They found 
a variety of sources or cores associated with the IRDCs. Most 
of these 'sources' are seen in absorption at 8.3 /mi, but occa- 
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sionally also in emission. ICarev et all l|2000b) found masses 
ranging from 20 M0 to 1000 M n and temper a tures from 
cold (10 K) to warm (40 - 50 K). iRedman etail §003), us- 
ing SCUBA, examined the G79.3+0.3 region and found young 
stellar objects to be associated with sub-mm observations at the 
edge of the IRDC complex. 

Because of their sheer sizes, high column densities and low 
temperatures, it has been suggested that IRDCs are the sites 
of the formation of massive stars. However, unlike low-mass 
star formation, little observational and theoretical information 
is available on the earliest stages of the formation of high-mass 
stars. Many high-mass star formation sites are known, but they 
are always in a relatively evolved state, where temperatures are 
> 30 K rather than 10 K. Therefore, if high-mass star formation 
can be linked to IRDCs, this would complement the evolution- 
ary picture. 

Low-mass cores have been studied in great detail in the 
sub-mm. Models of the spectral energy distribution allow 
to constrain the density and temperature structures of the 
cores to be compared with theoretical models. These observa- 
tions generall y suppor t theories of star formation in molecu- 
lar co res JShirlev et allkOOOllEvans et all 1200 ltlShirlev et all 
120021) . Similar techniques have been applied to high-mass pro- 
tostars in order to probe the physical structure of the en- 
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Fig. 1. InfraRed Dark Cloud in the direction of W51 as ob- 
served by MSX with SCUBA contours overlaid. The MSX 
image is stretched from 2.0 x 10~ 6 W rrT 2 sr _1 (black) to 
4.5 x 10" 6 W m 2 sr 1 (white). Overlaid are SCUBA emission 
contours at 4,6,8 and 10 times the 850 fim rms noise level of 
o"850 = 0.040 JybeairT 1 . The zero point of this image is at 
19'>21 m 44.7 s , 13°49'34.7" (J2000). 



velopes surrounding these sources ( van der Tak et al. I l2000t 
lHatchell & van der TakLEoollMueller et allkOO^ .The differ- 
ence between the observations described here and those listed 
above, is that all but two of the sources listed in the literature 
have strong mid- to far-IR emission which indica tes the sources 
are rel atively well evolved. The two sources of iMueller et alJ 
d2002h . which show no mid-IR emission, were only observed 
at 350 fim and therefore not modelled. 

There are many questions regarding IRDCs. Do they repre- 
sent the earliest stages of high-mass star formation? If so, how 
do they compare with Bok globules and low-mass star forma- 
tion sites? If not, what is their fate? In order to begin address- 
ing these questions, SCUBA 450 fim and 850 fim data have 
been obtained toward one IRDC complex seen in the direction 
of the W51 GMC. Furthermore, HCO + line observations were 
taken. This article presents and discusses the physical proper- 
ties (density and temperature structure and internal heating) of 
the complex through detailed modelling of both the sub-mm 
dust continuum emission and the HCO + line data. 

Section |3 discusses the observational data and the data re- 
duction. In Sect. the SCUBA and MSX observations are 
compared to estimate the optical depth. In Sect. [3] the contin- 
uum data is modelled using a spherically symmetric core model 
with the aim to clarify the density and temperature structures. 
These physical structures are used as input for the modelling 
of the HCO + line data (Sect.|4j, from which constraints on the 



velocity structure and the dumpiness are derived. Finally, an 
overall discussion is presented in Sect. ??. 

2. Observations 

2.1. MSX Data 

This investigatio n has mad e use of MSX Galactic Plane sur- 
vey data JPrice et allboOll) . The MSX radiometer, SPIRIT III, 
has surveyed the Galactic Plane in four mid-IR bands with 30 
times higher spatial resolution than IRAS. The MSX spectral 
bands are centred at 8.3, 12.1, 14.7 and 21.3 fim. The 8.3 fim 
spectral band is especially sensitive and covers the brightest of 
the Galactic PAH features between 6-9 fim. It is against this 
Galactic emission that IRDCs are seen. 

The IRDC considered in this article is a complex clearly 
visible in MSX images in the direction of the W51 GMC (Fig. 

The image is in equatorial coordinates and has been cen- 
tred on the peak extinction which is also at the centre of the 
S-shaped extinction valley. The greyscale is chosen such that 
black regions correspond to areas of lowest 8.3 fim emission, 
or highest extinction by the IRDC. Section l2~4l elaborates on 
the MSX extinction and on its relation to the SCUBA 850 fim 
emission, which is shown by the overlaid countours in Fig.^ 

2.2. JCMT Continuum Observations 

JCMT 1 observations of the W51-IRDC were made on October 
28, 2002 at 450 fim and 850 fim with the Sub-millimeter 
Common User Bolometer Array. The SCUBA observations had 
good weather conditions with an atmospheric optical depth at 
225 GHz of tcso ~ 0.06. To cover the extinction region seen 
in MSX, a 6' by 4' area was observed using the scan map 
mode. Two perpendicular chop angles were used with three 
chop throws of 30", 44", 68" a nd the image was rec onstructed 
using the Emerson II technique djenness et al.Ll20 00i). With the 
Emerson II image reconstruction technique, spatial structures 
much lo nger than the 68"chop throw are not reliably recon- 
structed djohnstone etall Eo03). This implies that a zero level 
correction is necessary and that this correction is local, i.e., it 
changes over the image on arcminute scales. 

The data were flat-fielded, extinction corrected, and cali- 
brated on_IJranusjisirig_fhe standard SCUBA reduction soft- 
ware djenness & Lightfooi [l998). During the observations, 
SCUBA was experiencing the transient noise phenomenon as- 
sociated with contamination caused by superfluid helium films. 
Over a few days, this film moved along the bolometer arrays, 
and in this case, moved from the 850 fim array to the 450 fim ar- 
ray. However, over the period of a few hours the noise was sta- 
ble enough to identify affected bolometers and to remove them 
from further processing. The removal of these bolometers had 
as main result the increase of the noise in the final image above 
the value expected for the assigned observing time. From the 

1 The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the parent organisations Particle Physics and 
Astronomy Research Council in the United Kingdom, the National 
Research Council of Canada and The Netherlands Organisation for 
Scientific Research. 
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Table 1. Fit to the SCUBA beam pattern as observed on 
Uranus. The total beam is approximated to be a superposition 
of Gaussians. The 450 fim beam strongly deviates from a sin- 
gle Gaussian, with more than 50% of the power falling into the 
error beams. The beam profiles are plotted in Fig. [5] 



Uranus at 850 fim Uranus at 450 fim 
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final images the noise was estimated to be 250 mJy beam 1 at 
450 fim and 40 mJy beam -1 at 850 fim. 

2.2.1. Calibration and Image Reconstruction 

Although the calibration should be good to 15% and 25% fo r 
the 850 W and 450 W filters, respectively (Jenne sset all EoOO). 
the observations were obtained early in the evening when the 
Flux Conversion Factors (FCF: the conversion from engineer- 
ing units of Volts to Jy beam -1 ) can change significantly dur- 
ing the time of the observations (Wouterloot, priv. comm.). For 
the observations presented in this article the FCFs are likely to 
have declined by about 10% and 20%, respectively, for 850 fim 
and 450 fim. Unfortunately, Uranus observations were not ob- 
tained before and after the IRDC observation, thus a potential 
systematic calibration error had to be taken into account. 

It is, however, possible to estimate the impact of this cali- 
bration change by re-performing the fit of Sect. l3.3l with the ar- 
tificial reduction of the 850 fim and 450 fim calibration factors 
by 10% and 20%, respectively. The main effect of this reduc- 
tion is to broaden the ^-topology, while the best fitting param- 
eters are not changed significantly. A varying FCF, although 
increasing the uncertainty in the calibration, does not alter the 
main findings of this work. 

2.2.2. SCUBA Beams 

Since both SCUBA 450 fim and 850 fim beam patterns are quite 
complex and the modelling was restricted to one dimension, az- 
imuthal averages of the Uranus observations were made (Fig. 

The actual beam pattern differed significantly from a sin- 
gle Gaussian, especially at 450 fim where 55% of the power 
was not in the main beam. The actual beam pattern, B(9), was 
approximated by a superposition of Gaussians, i.e., 



B(ff) = wxG x (6, <Ti) + w 2 G 2 (0, cr 2 ) + . . . , 



(1) 



where the Gaussians were normalised to account for the two 
dimensional nature of the beams and are given by 



G(6,cn) = 



1 



2tto-; 



exp 



"2cr? 



Fig. 2. Beam profiles at 850 fim (left) and 450 fim (right) as 
measured on Uranus on October 28, 2002. The beam profile is 
represented by a superposition of Gaussians, whose parameters 
are listed in Tabled The beam profiles are normalised to their 
central (6 = 0") value. 



(2) 



The relative weights w, of the Gaussians and their FWHMs 
(= 2.35 cr) are listed in Tabled Three Gaussians were used 
to characterise the 450 fim beam, while two Gaussians were 
sufficient for the 850 fim beam. These approximated beam pat- 
terns were used in the interpretation of the observed maps in 
terms of the modelling as described in Sect. 13.21 



2.3. JCMT Line Observations 

Observations of the W51-IRDC region were obtained in a va- 
riety of molecular lines. With JCMT, maps were obtained in 
HCO + in both /= 3^2 (267.6 GHz) and /= 4^3 (356.7 GHz) 
transitions on 15 September 2003. These maps were made on 
a rectangular grid spaced 5" in Right Ascension and 10" in 
declination. The maps are shown and discussed in Sect. |4] 
Moreover, single pointing observations of HCO + 3— » 2 (dated 
December 5, 2001) and H 13 CO + 3^2 (dated June 15 and July 
22, 2004) are available at two positions. Main beam widths are 
20" and 14" with efficiencies of rj = 0.69 and rf = 0.61 for the 
3— >2 and 4— >3 transitions, respectively, and the velocity reso- 
lutions are 0.35 km s -1 and 0.53 km s . Furthermore, similar 
to the 850 fim continuum case, it was assumed that the beams 
could be approximated by two Gaussian components where the 
power of the error beam is given by 1 - r\ and its FWHM is de- 
rived from the SCUBA's 850 fim beam by assuming a scaling 
proportional to A. 

All spectra contained low frequency standing waves and 
a linear baseline subtraction was performed in two windows 
on opposite sides of the line, covering 15-25 km s -1 and 
41-51 kms -1 . The IRDC emission is centered at v'ljjr - 
34.0 km s -1 with a spatial variation of ~ 0.5 km s -1 over the 
complex. Using the description of the galactic rotation curve 
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W51 IRDC detailed view 
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Fig. 3. This scatter plot shows the anti-correlation between 
8.3 fim emission and 850 fim emission. To model the trend, 
a median of the MSX emission within bins of 0.02 Jy beam -1 
was measured. The squares and associated error bars indicate 
the median values and upper and lower quartiles. The best fit- 
ting model to the median data is indicated by the solid line. 



PEAK VALUES 



Table 2. 450 fim and 850 fim positions, peak intensities and 
8.3 fim optical depth of PI, P2 and the extinction peak. 



bv lWouterloot & Brand! J 1989) this can be translated into a dis- 
tance of 2.7 kpc. 2 

2.4. Mid- Infrared Extinction 

With images of the region in both the mid-IR and the sub- 
mm it is possible to estimate the 8.3 fim opacity within the 
dark cloud complex, together with the background and fore- 
ground e mission of the G alactic Plane. Following the approach 
of Ijohnstone et alJ d2003l) . the measured anti-correlation be- 
tween 8.3 fim and 850 fim is interpreted as a column density 
measure in terms of a dust model with constant opacity ratio. 
Here it is assumed that the foreground and background intensi- 
ties at 8.3 fim, 7f g and 4 g , are constant over the map, such that 
the variation of the measured 8.3 fim intensity can be attributed 
to the IRDC opacity. When a constant temperature and optical 
thinness at sub-mm wavelengths are assumed for the IRDC, 



2 The distance assignment from the radial velocity suffers from a 
basic ambiguity. However, the far distance solution of 8.6 kpc can be 
excluded based on the appearance of the IRDC as dark foreground 
material in front of the bright MIR emission from the inner Galaxy. 
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Fig. 4. (top) Detailed view of the emission peaks, showing PI, 
P2 and the extinction peak at 850 fim. This figure is a close- 
up of Fig. ^ with the same contour levels, (bottom) 450 fim 
contour levels drawn at 4,6,8 and 10 times the mean rms-noise 
with cr 450 = 0.235 Jy beam -1 . 



this opacity, Ts.3 scales linearly with the 850 fim emis- 
sion, T8.3 ^ m = kl^o fan. The measured anti-correlation between 
^8.3 iim and 7850 fxm can then be used to fit the values of 7f g and 
7b e with which the 8.3 fim opacity can be computed. Because 
this fit assumes isothermal and constant dust properties, all po- 
sitions inside two apertures with radii of, respectively, 27" and 
18" centered on the two brightest 850 fim peaks had to be ex- 
cluded from the derivation of the foreground and background 
intensities as they clearly belong to sub-mm sources. The re- 
sulting scatter plot is shown in Fig. [21 

To measure the trend more distinctly, the data were binned 
in 14 bins of 0.02 JybearrT 1 . The solid squares were calcu- 
lated as the median 8.3 fim emission within each bin. The error 
estimates are the upper and lower quartiles around this median. 

The best fit to the trend gave a foreground emission of 
7f g = (1.6 ± 0.5) x 10~ 6 W m sr , a background emission of 
7b g = (2.3 ± 0.4) x 10 6 W irT 2 sr~' and a proportionality ratio 
of k — 5.1 + 2.7. This is shown in Fig.[3]as the solid line. 

Assuming a typical dust temperature of 15 K for the 
IRDC, the £-ratio translates in an mid to far-IR opacity ratio 
of ~ 870. This ratio is consistent with the va lue of ~ 640 
found for another IRDC ijohnstone et all 120031) and is close 
to th e value of ~ 500 predicted for the ice mantle dust grains 
bv lOssenk onf & Henning ( 1994). In contrast, the diffuse ISM 
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grains of iLi & Draind J200ll) have a considerably larger ratio 
of - 1700. 

With the fits to the foreground and background emission, 
a rough estimate for the 8.3 fim, MSX-beam averaged, optical 
depth through the IRDC can be made, i.e., 

I -I f . 

(rg.3fim) = -In— : . (3) 

Where / is the measured intensity from the MSX-map and I\, g 
and 7f g are the estimates from the fit above. The 'MSX extinc- 
tion' therefore anti-correlates with the measured 8.3 fim inten- 
sity. 



3. Continuum Modelling 

3.1. Radial Profiles 

The 450 fim and 850 fim data show three distinct emission 
peaks (see Fig.0}. The emission peaks are labelled as PI, P2 
and EP and the coordinates are listed in Table[2] Table|2]also in- 
dicates the peak brightness and the estimated extinction as ob- 
tained from Eq. Q. The peak of MSX extinction (position EP) 
clearly corresponds to an 850 fim emission peak. At this posi- 
tion 450 fim also peaks, but it is relatively weak. In contrast, 
the other sub-mm emission peaks are strong at both 850 /im 
and 450 /urn. Note that PI lies at the edge of the MSX extinc- 
tion structure, while P2 and EP are well within the complex. 

In order to model the emission, one-dimensional radial pro- 
files were c onstructed from the sub-mm data. Following the 
approach of Muel ler et alJ J2002I) . the sub-mm images are az- 
imuthally averaged over relatively clean regions. The extended 
IRDC emission was avoided by choosing azimuthal angles for 
which a radial cut does not pass over other parts of the cloud. 

Figures [5] and [6] show azimuthal averages of the intensity 
around the three peak positions. Here a zero-level radius of 60" 
has been introduced such that all profiles go to zero intensity at 
this radius. In this way low density material at larger scales is 
ignored, because no reliable information is available for these 
scales due to the limited chop angle. It is clear that the choice 
of the zero-level radius introduces systematic uncertainties in 
the parameters derived from the model. As such the derived 
radius of the core, R, has to be smaller than the zero-level radius 
and the derived core mass, M(R), is limited to the mass well 
within this zero-level radius. From tests with varying radii we 
found, however, that the phyiscal structure, i.e., enclosed mass 
as function of radius, M(f), and temperature structure, T(r), 
are quite insensitive to the applied zero-level radius. The 60" 
zero-level cut-off applied here is the most natural choice for 
the boundary between the core and extended material in the 
IRDC. At these radii the intensity profiles flatten out, indicating 
that the emission originates here from diffuse, extended regions 
rather than the core. One has to keep in mind, however, that 
the derived total mass of the cores is essentially constrained by 
this definition of the zero-level radius, as there is no sharp outer 
boundary of the cores. 



3.2. Core Model 

Each position, PI, P2 and the extinction peak was modelled as 
a spherical core with a radial density distribution falling off 
as a power-law. The cores are heated by the interstellar ra- 
diation field (ISRF) and a possible internal radiati on source. 
The on e dimensional diffusion code CSDUST3 of Egan et alJ 
dl988l) was used to calculate the continuum emission from each 
of the cores. Together with a density power-law exponent, the 
code requires a dust grain model that provides the opacities, 
the optical depth through the core at a specific wavelength, the 
intensity of the ISRF and the luminosity of the internal source. 

CSDUST3 allows for a wide range of adjustable parame- 
ters. For reasons of efficiency many of these parameters have 
been fixed throughout the modelling: 

- Because the data show all three cores as roughly circular, 
only spherical geometries were considered. Although on 
a larger scale the cores are part of an elongated extended 
structure, no attempt was made to incorporate this geome- 
try in the modelling; 

- The dust model of Ossenkorjf & Hennin 3 (GUI model 5; 
OH5) was used. This reflects a high density environment 
where intermediate ice-mantles might have been formed. 
These opacities are consistent with the extinction analysis 
of Sect. I2.4I OH5 opacities are given in units of cm 2 per 
gram refractory material and the dust-tq-hyd rogen ratio of 
1.5 x 10~ 26 g H-atom" 1 dDraine & Leell 19841) was applied 
to convert the opacity to units of cm 2 per hydrogen atom. 

- To avoid infinite densities an inner cavity radius r c was 
introduced. The radius of the cavity was chosen to be 
r c = 0.001 pc = 206 AU which is small compared to the 
outer radius of the core (~ 1 pc). Hatchell & van de r Tak 
J2003I) concluded for similar modelling that their density 
distribution and resulting intensity profiles were insensitive 
to the radius of the inner cavity. 

- The wavelength grid consisted of 60 points, distributed be- 
tween A = 1 fim and 1 300 fim. The spatial grid consisted of 
100 points and was chosen to have a density of grid points 
proportional to r~ 2 . 

- The model makes use of a 'reference frequency' to which 
the opacity is scaled. 450 fim was chos en as the reference. 

- A description of the ISRF is given bv lMathis et ail ril983h 
and their ISRF at a Galactocentric distance of 6 kpc was 
used. Since the wavelength grid only starts at A = 1 fim, the 
UV and part of the near-IR components have been left out. 
This is justified since the cores lie well embedded within 
the molecular cloud. 

- In cases where internal heating was included, an effective 
temperature, T e ff, of this 'heating' was fixed at 6 500 K 
for all models. A higher value of r e ff had no effect on 
the temperature structure, since the high energy photons 
are quickly re-radiated. The emergent profiles are only af- 
fected when the T e g of the inner heating becomes very low 
(r eff < 30 K). 

With these specifications four parameters remain free to op- 
timise: 2) the exponent of the density distribution p, ii) the op- 
tical depth at the reference frequency r ro( j, Hi) the size of the 
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Fig. 5. 850 fjm radial profiles of PI (left panel), P2 (centre) and the extinction peak (EP, right). For each panel three features 
are shown. First the radially averaged data points, binned at every half main beamwidth, where the radial averaging involves 
directions not contaminated by excess emission from, e.g., neighbouring cores. These points were normalised to their r = 0" 
values. Second, the dashed line shows the beam profile as obtained from the Uranus observations, which was also normalised to 
its central value. Finally, the solid line gives the best model fit to the data (see Table[3). 
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Fig. 6. Same as Fig. |5] but for the 450 pm data. 
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cloud, i.e., the outer radius, R, and z'v) the luminosity of the 
heating source at its centre, L. 

The first three combine to give the mass of the core 

M = 4™ H £^^^£* 2 (p*l,3), (4) 
2 3 - p cr rod A ^ _ i 

where is the mass of one hydrogen atom, p/2 = 1.18 
the mean molecular weight per H-atom, <x rod = 1.01 x 
10~ 25 cm 2 H-atorrT 1 the opacity at 450 pm and A DR = R/r c is 
the dynamic range. Because the mass of the core is as a phys- 
ical quantity of more interest than the reference optical depth, 
M was chosen as an input parameter and T ro d was adjusted ac- 
cording to Eq. ©. 

The output intensities of the model as functions of impact 
parameter were convolved with the approximated beam pattern 
from Tabled For each set of parameters a^- 2 goodness-of-fit 
was calculated 




with / the observed intensity, M the convolved intensity of the 
model and cr the error in the observation. The respective con- 
tributions to the^ 2 denote, respectively: 



- The 450 pm and 850 pm profiles, which contribute the data 
points of the profiles within the cut-off radius of 60", sepa- 
rated by the main beams' half width. 



- The IRAS non-detection at 100 pm. M\qo is the 
model intensity at 100 pm convolved with the IRAS 
beam (approximated to be spherical with FWHM = 220") 
and o"ioo is the rms-noise of IRAS, which is estimated at 
o- 10 o = 30MJy sr -1 . 



3.3. Results and Analysis 

To find the x 2 minimum in the 4-dimensional parameter space 
the problem was reformulated as an inverse problem for which 
soluti on (or retrieval) techniques are available (e.g.. lRodgerst 
2000). An iterative scheme was applied and the absolute min- 
imum of x 1 was usually found in a few iterations. Moreover, 
X 1 contour plots were made in which two parameters were var- 
ied. In this way a good visual impression on the^ 2 topology is 
obtained. 

Table|3gives for every profile the^- 2 minimum and its loca- 
tion in the (p, At, R, L)-plane. The errors in the parameters give 
the maximum range for which x 2 increases by 1, with the re- 
maining three parameters free to find their optimum values, i.e., 
by minimising x 2 - This Ax 2 b ehaves also as a redist r ibutio n 
with one degree of freedom dBevington & Robinson! Il992h . 
The Ax 2 = 1 surface gives therefore the lcr confidence level 
of the parameters. Finally, the (MSX beam averaged) model 
optical depth at 8.3 pm is calculated to provide a comparison 
with the MSX results. 
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Table 3. Results of the profiles that have been fitted. In Cols. 
(2-5) the optimum values of p, M, R and L are given for which 
X 2 (Col. 6) is at its minimum. Indicated in the sub- and su- 
perscript of the parameters is the range of this parameter for 
which the^- 2 increases by 1, indicating the lcr-limits. Finally, 
the modelled beam-averaged rg.3 ^ m (Col. 7) is given. 



X 2 Contourplot 




log Core Moss [M s ] 

Fig. 7. x 2 contour plot of the PI model with fixed density ex- 
ponent (p = 2.25) and diameter (75")- The solid lines indicate 
the value of^ 2 in the M - L plane. The elongated shape of the 
X 2 contours illustrate a basic mass-luminosity degeneracy. The 
star indicates the x 1 minimum which is close to the absolute 
minimum given in Tabled 

Examples of x 2 contour maps are given in Figs. and [S] 
In Fig. [7] x 2 contours are plotted as functions of mass and lu- 
minosity at a fixed size and density exponent, while in Fig. [S] 
density exponent and luminosity are the varied parameters. 

We can derive some general conclusions from the^- 2 maps: 

- The IRAS detection limit, estimated at 3cr =s 90 MJy sr -1 , 
corresponds roughly to luminosities of L — 500 - 1 000 Lq, 
only weakly dependent on mass or density exponent. 
With luminosities approaching 10 3 Lq, the contribution of 
100 /im to the total x 2 becomes significant. This sets upper 
limits to the luminosities of the embedded sources. 

- For the PI and P2 models a mass-luminosity degeneracy is 
observed. Figure gives an example, where both a At = 
130 Mq, L = 150 L -model and a M = 70 M Q , L = 
600 LQ-model yield x 2 = 15. Core mass and radius are 
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X 2 Contourplot 




Density Exponent 



Fig. 8. Another contour plot of the PI model, now at a fixed 
mass of 90 Mq and a diameter of 75". The range in luminosity 
for which good fits (x 2 < 15) can be obtained spans almost 
one order of magnitude, while the density exponent can vary 
between 2.0 and 2.5. 

the quantities which are best constrained by the continuum 
data. 

- Similar to the M - L degeneracy a p - L degeneracy can 
be observed (see Fig.|S|l. A reduction of the luminosity by 
a factor of three can be compensated by the increase of the 
exponent by 0.2. 

- For the extinction peak all x 2 levels are lower than for the 
other two peaks. This is, however, only due to the large un- 
certainties that have been attributed to the 450 pm profile. 
It is the 850 pm profile that dominates the^ 2 -fit. The error 
bars on the parameters are actually comparable between the 
cores. 

- For PI and P2, on the other hand, the main contribution to 
X 2 originates from the 450 pm data. 

The density and temperature structures that result from the best 
fit parameters are shown in Figs. 1^1 and 1 1 01 For the extinction 
peak temperatures have a minimum of about 9 K within the en- 
velope of the core. At the outer boundaries the temperature is 
set by the intensity of the ISRF and at the inner boundary by the 
19 Lq luminous source. The other two peaks (Fig. HOt. though 
with higher absolute temperature levels, follow the same qual- 
itative structure. However, since the lower luminosity limit of 
the extinction peak is only 1 .7 Lq, a source-free model is not 
completely discounted and its density and temperature struc- 
ture is also given in Fig. [9] Here, we have only a marginal indi- 
cation of an internal heating source from the continuum obser- 
vations. 

4. HCO + Line Modelling 

4.1. Lines 

The results of the continuum model above provide constraints 
on the density and temperature structures of the three cores. 
Modelling of HCO + lines gives information on the velocity 
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Fig. 11. Contour maps of HCO + 3— >2. (left) Integrated intensity in a window of v = 34.0 + 2 km s _1 . Contour levels range from 
1.0 K km s _1 to 3.0 K km s _1 in intervals of 0.5 K km s _1 . The orientation is the same as in the continuum images. The positions 
of the peaks and the JCMT beamwidth are also indicated, (right) The intensity of the highest channel in the v = 34.0 + 2 km s _1 
window. Here the contours start at 0.70 K with intervals of 0.2 K. 
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Fig. 9. The solid lines give the temperature (grey) and density 
(black) structure of the minimum x 2 model of Table [3] for the 
extinction peak. It shows that temperatures drop below 10 K 
within the envelope, but rise toward the boundaries through the 
heating of the ISRF and the 19 Lq source, respectively. The 
dashed line represents a model without inner heating (L = 0). 
Here, temperatures drop below 5 K and densities are higher 
to compensate for the lower temperature. In reality cosmic ray 
heating would have stabilised the temperature, but the differ- 
ence would not be noteable from the observations as it is con- 
fined to small scales. The x 1 of this model is 3.9 and its core 
mass equals M = 170 Mq. 



structure and, using a turbulent-clumping model, puts con- 
straints on density inhomogeneities. The advantage of using 
HCO + is that due to its high Einstein A coefficient and the 
high densities i n the cores it radiates easily. Moreover, theo- 
retical studies o f lRawlings et all J 19921) and observations of the 
Orion molecular cloud ( Savv a et al. 1 120031) show that HCO + is 
unaffe cted by the depletion of CO. Recent models of lLee etaD 
ll2004l) . on the other hand, claim the opposite. If HCO + is in- 
deed depleted over a wide region of the core, the total abun- 
dance will be higher than derived here. However, rather than 
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Fig. 10. The temperature (grey) and density (black) structure 
for the minimum^ 2 model of PI (solid lines) and P2 (dashed 
lines), respectively. Parameters are given in Table|3] 



deriving abundances, the main goal of the line model is to test 
whether the physical temperature and density structure as de- 
rived from the continuum can be confirmed, rejected or fur- 
ther constrained by the line fits. Moreover, we try to extract 
dynamical information from the lines. A disadvantage of using 
HCO + is that, contrary to the sub-mm continuum, the lines are 
optically thick (the models showed typical line-centre optical 
depths of t st 10), and the models are especially sensitive to 
the core's boundary layer. 

The line maps of the HCO + 3— >2 transition are shown in 
Fig. 11. The emission traces the general structure of the fila- 
ment. The three cores are best visible in the maximum intensity, 
but the overall continuum emission is only matched if the in- 
formation from the maximum intensity map and the integrated 
intensity map is combined. From the comparison of the two 
maps it is seen that the lines are narrow at the extinction peak 
and quite broad at P2. These differences in the lines shapes 
are exactly reproduced in the following fits of the line profiles. 
Perhaps the most notable difference between the two maps is 
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" Data as extracted from map or pointed observations. 
Table 4. An overview of the available lines for every peak. The 
data have been obtained during several JCMT observing runs 
(Sect. EH. 



the region east of P2, which is bright in the integrated map. In 
this region the lines are quite broad. As this broad line region 
is not associated with a sub-mm core we have not further anal- 
ysed its properties here, but will be the topic for future studies. 
It must be emphasised, however, that due to the noisy nature 
of the spectra, the rather large optical depths and the low res- 
olution, not much information can be obtained from the line 
maps. Thus no further quantitative conclusions are drawn from 
the maps apart. 

An overview of the available lines (including pointed obser- 
vations) is given in Table 4 and the observed lines are shown 
in Fig. 13. The S /N ratios of the lines reflect the heterogene- 
ity of the observations (see Sect. 12. 3> . For the extinction peak 
the HCO + 3— >2 pointing has a high S IN, while the isotope is 
barely visible. The HCO + 3— >2 transition of PI shows a nar- 
row sub-feature near v +2.0 km s _1 , which also appears in 
the isotope. At HCO + 4— >3 it is, however, absent. The shape of 
the sub-feature is suspect as it might be an instrumental artifact 
or blending with another cold clump. As the nature of the sub- 
feature is not understood it is assumed at least unrelated to the 
core and the corresponding channels are masked. However, the 
assumption that it is merely noise is also tested. In that case the 
derived parameters change very little relative to the 'masked- 
channels model' with the variation well below the error bars 
given in Table 5. 

The model, explained in more detail below, is structured 
similar to the continuum model. For a given set of parameters it 
computes a line spectrum, which, after the convolution with the 
appropriate beam structure (see Sect. 12. 3> . is compared to the 
observed spectral lines. In an iterative fit procedure the model 
parameters providing the best match to the lines in terms of a 
X 2 criterion are obtained, x 2 is defined as 



X 



y - M u 



i-j 



(6) 



where /y is the intensity in channel i of transition j, M the mod- 
elled intensity and <x the rms-noise in the spectrum. One run of 
the model includes both HCO + transitions and in the case of 
PI an d the extinction peak al so the H 13 CO + isotope. A value 
of 70 llWilson& Rood! 1 1994 is used for the HCO + /H 13 CO + 



SIMLINE DENSITY, TEMPERATURE STRUCTURE 
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Fig. 12. Schematic picture of the density, n(r), and temperature, 
T(r), structure. For SIMLINE these must be given as power- 
laws and the temperature structures of Figs.|5]and^|is mim- 
icked using two shells (except for the EP-c model, where one 
shell suffices). The first shell ranges from r\ to ri. Here r\ rather 
than r c is the inner radius and r2 is the position where the tem- 
perature is at its minimum. The second shell ranges from r2 
to R, the core radius. The values given here are those of the 
clumps and are indicative of PI. 



abundance ratio. For each peak between 28 (EP) and 32 (PI) 
channels are then simultaneously fitted. 



4.2. Simline 

The calculation of the model profiles was p erformed by the 
spheri cally symmetric SIMLINE program of lOssenkopf et all 
(2001). This code is unique since it takes into account effects 
of turbulent clumping, characterised by a turbulent velocity 
distribution and the turbulence correlation length /corr- The 
clumping, which can result in a significant reduction of the 
effe ctive opacity, is treated in a statis t ical a pproach following 
iMartin et all dl984 . lOssenkopf et alJ d200 lb have shown that 
this results in a reduction of the strong self-absorption effects 
from which all (homogeneous) microturbulent models suffer. 
The code uses a shell structure with parameters that vary as 
power-laws within the shells and that can be arbitrary from 
shell to shell. Following the results of the continuum, two shells 
were used to reflect the temperature structures of Figs. |5] and 
1101 as close as possible (see Fig. 12). Since the (clump) den- 
sities are high, the gas and dust temperature are expected to 
be well coupled. At the core radius R, however, the gas tem- 
perature was slight ly augmented to 20 K to acc ount for PDR- 
heating of the gas (Sternberg & Dalgarno, 1989). Assuming a 
shielding by the surrounding wide-spread low-density material 
of A v ~ 0.1 - 0.5 the PDR-heating effects produce only a few 
degrees difference between gas and dust temperature, so that 
the assumption of 20 K dust temperature should be accurate 
within about 1 K. The density structure, determined by p, is 
directly adopted from the continuum (Table [3j. Because it was 
found that the detailed structure in the centre of the core at radii 
r < r\ — 0.01 pc had no significant influence on the resulting 
line profiles, we have omitted this part in the radiative transfer 
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Fig. 13. HCO + observations (histogram) together with the best fits of Table 5 (solid lines). The three rows correspond to the 
transitions (from top to bottom: HCO + 3— »2, HCO + 4— >3, H 13 CO + 3— >2), while the horizontal direction indicates the model (PI, 
P2, EP, from left to right). No H 13 CO + data is available at P2. Furthermore, each panel shows the bounds of the fitted spectral 
window (dashed vertical lines), the rms-noise level (upper left) and in case of PI the masked channels (crosses). Velocities are 
shown relative to vlsr = 34.0 km s _1 . For comparison the EP-plots also show the best fit L — model with fixed velocity 
exponent, v exp = 0.5 (dashed line). 



computations. This increased the speed of convergence consid- with the continuum. For each of these models, four input pa- 
erably. rameters were varied: 



Four types of models were considered: For each of the 
peaks the model that corresponds to the parameters of the best 
fit to the continuum (denoted PI, P2 and EP) and, in the case 
of the extinction peak, the additional best fit source-free model 
(denoted EP-c). In this way the line model is self-consistent 



- The turbulent velocity structure (2 parameters: central 
value and radial exponent). Clouds are often o b served 
to have exponents of 0.5 (e.g., iMiesch et all [l999; 
Ossenkorjf & Mac Low, 120021) . while negative exponents 
indicate that the turbulence is driven from the interior. 
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- The combined mass. The program simulates a clumpy 
medium which comprises two components: The (by defi- 
nition identical) clumps, which fill a volume fraction / and 
the interclump region, which is treated as void. The density 
of the clumps is then the average density multiplied with 
the inverse filling factor 1//. Consequently, the combined 
mass, yVt com = Al//, is representative of the conditions 
inside the clumps. By fitting this parameter for the effec- 
tive density 3 of the clumps, and comparing the result with 
the mass derived from the continuum the filling factor / is 
obtained. 

- The mass of tracers. In the statistical treatment of the tur- 
bulent clumping, it is mathematically equivalent for the 
radiative transfer problem to treat the clumpy medium 
as a homogeneous structure with the properties of the 
clumps, though with a molecular abundance decreased by 
/. Therefore, the output is only sensitive to the combined 
abundance, X com = Xf. As the filling factor is already fit- 
ted via the combined mass, the molecular abundance of 
HCO + can be derived directly by joining the combined 
mass and combined abundance into the mass of tracers, 
MX = yvi com X com , which we use as the fourth free param- 
eter. 

Varying these four parameters is sufficient for the model to 
fit the data. Therefore, the radial gradient for the velocity struc- 
ture v ex P is assumed to be the same for both shells, and X com 
is assumed to be constant over the whole core. Other parame- 
ters that are fixed are the systematic velocities (e.g., the spectra 
do not show evidence for infall) and the turbulence correlation 
length, which is fixed at the typical value of /corr = 0.01 pc. 
For a more t horough discussion on these radiative transfer pa- 
rameters see lOssenkopf et aO d200ll) . 

4.3. Results 

Results are presented in Table 5, where the best fit parameters 
are given. The error limits denote the 1 cr boundary of x 1 , 
calculated in the same way as the continuum. The 'new' results 
are then y turr \ y ex P, X and /, though X and / are not not 
directly fitted but derived from the fit parameters using the 
mass values from the continuum fits. The error limits on / 
and X are calculated by appropriately dividing the upper and 
lower error limit, e.g., X up = (XM) up /M low , where M low is 
the lower 1 cr mass-limit of the continuum. 

All fits are satisfactory, as can be seen in Fig. 13, which 
demonstrates the agreement between fitted and measured pro- 
files. The agreement can be quantified using reduced x 2 val- 
ues x\ - X 2 /(number of frequency bins - number of free 
parameters) (see Table 5). Except for the cold extinction peak 
model they fall below one and even for this model the value is 
only slightly above one. 4 Thus the line fit shows again that the 
model with the internal heating fits the data better, though the 

3 _yvf om is merely a scale function for the density, with the radial 
gradient fixed by the continuum model. 

4 If the 'subfeature' in the PI lines is not masked, the xl raises to 

1.4. This results only in slight changes of the derived parameters. 



cold model still gives reasonable fits. It is remarkable that the 
errors on the parameters are similar for PI and the extinction 
peak in spite of the higher S/N for PI. This can be explained 
from the equivalent total amount of information for both cores. 
As the absolute height of the error bars is similar, the relatively 
weak emission in the isotope for the extinction peak constrains 
the model in the same way as the stronger emission for PI. 
In contrast, the error bars of the parameters derived for P2 are 
considerably larger due to the lack of H 13 CO + data. When con- 
sidering the general topology of the x 2 surface in the consid- 
ered four-dimensional parameter space in all cases a clear and 
prominent global minimum was found, similar to the contin- 
uum case. By considering many different initial guesses and 
comparing the results it can be guaranteed that no other local 
minima were accidentally hit. 

No explicit test was made how the uncertainties in the fit 
results from the continuum modelling affect the outcome of 
the line fitting. A general estimate for this sensitivity can be 
obtained, however, from the fit results for the extinction peak. 
Here, we compared a hot and a cold model for the same core. It 
is found that, although the physical structures of the two mod- 
els are notably different (see Fig. [5}, all derived gas-specific pa- 
rameters are very similar. We therefore anticipate that the errors 
in the dust model would hardly propagate into larger noticeable 
errors of the corresponding HCO + models. 

Because for the sourceless EP-c model the line fit is not 
expected to yield negative values for v exp , we have explicitly 
tested its influence on the EP-c model by fixing the exponent 
to and 0.5. This caused the x 2 of the best fitted model to 
increase to 36.1 and 39.1, respectively. The increase is due to 
a slight narrowing of, especially the 4— >3 line. Thus, due to 
the negative v exp , the turbulent velocities at the core's centre 
are large and allow photons to escape through the line-wings, 
thereby broadening the emergent line-profile. The increase of 
X 2 by more than one can be taken as a statistical confirmation 
that it is valid to fit the v exp parameter and that the error limits 
given in Table 5 are significant. Fig. 13 shows the difference 
between the two self-consistent models for the extinction peak: 
The 'hot' model where v exp is allowed to find its optimum and 
the 'cold', L = 0, model with v exp fixed at 0.5. The^ 2 difference 
of * 12, arises mainly due to the 4 — > 3 profile, which the 'hot' 
model fits much better. 

5. Discussion 

Although the pure ^-values do not exclude the cold model, 
it can be excluded based on considerations with respect to the 
dynamical structure of the core. All model fits provide nega- 
tive values for v ex P, i.e., a decrease of turbulent motions out- 
wards of the core. This clearly suggests that the turbulence is 
driven by star formation activity in the interior of the core. The 
same source that heats the dust and gas in the centre, is driving 
turbulent motions which are traced by the HCO + line profiles. 
From the combination of the continuum flux profiles and the 
line ratios, which give both only a relatively weak indication 
of a heating source in the core, with the line profiles indicat- 
ing a driving source, it becomes clear that also the extinction 
peak hosts an active young stellar or protostellar object. The 
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Table 5. Results of the HCO + line fitting routine. Column (1) shows the four models, the first three correspond to the best fit 
results of the continuum, while 'EP-c' corresponds to the best fitting L = model for the extinction peak. Cols. (2-5) give the 
optimum values of the free parameters at which x 2 (Col. 6) is minimised. The normalised^ 2 value are given in Col. (7). The 
error limits have the same meaning as in the case of the continuum. Comparing Cols. (4) and (5) with the mass obtained from the 
continuum model then gives the filling factor, /, and the HCO + abundance, X. 



cold model can thus be excluded. For the two sub-mm cores 
this is much more obvious. All parameters point consistently 
to luminous sources in both of them. 

A significant difference between the fit parameters of the 
cores is also seen in the molecular abundance. For the PI 
and EP models HCO + abundances of around 1.5 x 1CT 9 
are found, in good agreement with the results presented by 
Helmich & van Dishoeckl ( H997i) for warm cores within the W3 
complex. Both EP models give about the same value indicating 
that the derived molecular content does not discriminate be- 
tween the changes in temperature and density structures of the 
two models. P2, on the other hand, shows much higher values 
for the abundance. Even the lower limit of 3.1 x 10~ 9 falls at 
about twice the average value of the other cores. The reason, 
nevertheless, for the high fitted abundance can be clearly seen 
in the line profiles. P2 shows by far the strongest signatures of 
self-absorption in the profiles, asking for high optical depth in 
the HCO + line. In contrast one would expect a similar abun- 
dance for the three cores, since they resulted from the collapse 
of the same molecular cloud. A possible explanation for this 
contradiction is the existence of extended cool foreground ma- 
terial which contributes to the self-absorption profile, but is not 
detected in the continuum observations due to the limited chop 
angle which asks for the artificial zero-levelling at large scales. 
It is also intuitively clear that this effect is strongest for P2 due 
to its central position in the IRDC. Alternatively, a more com- 
plex temperature structure or a systematic change in the inter- 
nal clumping could also result in an enhanced self-absorption 
effect. One has to take into account, however, that at least a 
considerable fraction of the deviation may be due to the large 
uncertainty in the determination of the abundance for P2. These 
errors would have been severely constrained if H 13 CO + data, 
being a good probe for the column density and thus for XM 
due to its low optical depth, were available on P2. Additional 
observations, both on large scales and with higher resolution, 
are required to solve this puzzle definitively. 

The values of 0.1 - 0.2 found for the filling factors in- 
dicate a slightly clumped medium , in agreement with val- 
ues for other star-forming sites (e.g., Stutzki & Guesten, 1990; 
lOssenkopf et all EoOl ). Often even lower values of / can be 



found (e.g.. lHogerheiide et all fl995). though in these models 
the interclump medium is not treated as void. Apart from the 
dumpiness no conclusions can be drawn on material at densi- 
ties n > 10 6 cm 4 , well above the the critical density of the ob- 
served lines. The model is, however, sensitive to densities much 
lower than n cl -i t . The line profiles provide a strong constraint on 
the magnitude of self-absorption, thus allowing to derive the 
structure of the core's outer layers. As a result the model fit in- 
fers from the relatively weak self-absorption that the material 
in the outer layers must be clumped to remain sufficiently ex- 
cited. For PI, the filling factor is lower by a factor of two and it 
is tempting to link this low value to the relative transparency of 
PI at 8.3 fim in the MSX data (see Table|2j>- In the macroturbu- 
lent limit a lowering of / by a factor of two would also decrease 
the effective optical depth by a factor of two. This is, however, 
an extreme and the translation of the stronger clumping into a 
higher 8.3 /urn transparency is probably less efficient. Thus the 
clumping alone cannot explain the difference between the op- 
tical depths in Tables |2]and[3] Note that the error bars given in 
Table|2]are systematic, i.e., the uncertainties in the background 
and foreground emission are the same for all peaks, making the 
average extinction in P2 and EP at least twice the extinction 
measured for PI, probably even higher. However, it is prob- 
lematic to compare the findings of Tables[2]and[3] The models 
deal with the cores only and do not account for the surrounding 
cloud and the continuum models do not take into account the 
internal clumping at all. In summary, the models have clearly 
shown that the medium must be clumpy on an unresolved scale; 
Table 5 shows that the line profiles cannot be explained with a 
homogeneous medium. The stronger turbulent clumping of PI 
gives a partial explanation why this core is less prominent in 
the MSX data compared to P2 and EP, but it does not explain 
the full magnitude of the difference. 

The key question that triggered this research was whether 
star formation is taking place within the IRDC-cores. This is 
obvious for PI and P2, but for the extinction peak the contin- 
uum model allows the sub-mm emission to be explained just 
by density enhancement. The error bars on the 450 yum profile 
of Fig. |6] are simply too large to break the p - L degeneracy 
(see Fig.|8}: The presence of a possible source can be 'hidden' 
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by the model by reverting to a steep density profile. However, 
the velocity structure, as derived from the HCO + line results, 
clearly indicates a 'hot' model. We have shown, by combining 
the continuum and line models, that the extinction peak could 
be one of the earliest stages of star formation observed so far. 

The results for the emission peaks, on the other hand, point 
convincingly to the presence of bright sources, although their 
luminosities are poorly constrained and could be easily off by 
a factor of 2-3. The model does not constrain the nature of the 
sources, though with luminosities of hundreds of solar units, it 
is hard to think of anything but stellar. Pre main-sequence evo- 
lution tracks of lPalla & Stahlerldl999l) show 300 L Q luminous 
sources to correspond to masses in the range of M 3 - 5 Mq, 
only slightly dependent on the evolutionary stage. This cor- 
responds to intermediate-mass star formation and the spectral 
type would be late B. However, since star formation in gen- 
eral is not an isolated process, multiple stars might have been 
formed, which could change this analysis, although it is the 
highest mass star that produces the bulk of the luminosity. 

6. Conclusions 

The main conclusions of this study are: 

- Through self-consistent modelling of dust emission the 
presence of warm material embedded in a surrounding en- 
velope is identified for the cores PI and P2. While the na- 
ture of this heating is still unresolved, the energy involved 
should be at least a few hundreds of solar luminosities. For 
the extinction peak the presence of a source is not inferred, 
although the continuum model does prefer it. 

- The molecular line model of HCO + supports the physical 
structure as obtained from the continuum. The HCO + data 
reveal the clumpy nature and also provide a strong indica- 
tion for the existence of an energetic source in the extinc- 
tion peak. Only from the combination of line and contin- 
uum observations can the presence of a star, albeit in a very 
early evolutionary state, be inferred. 

- The derived luminosities of ~ 300 L© for PI and P2 are 
indicative of intermediate-mass star formation (M 3 - 
5 M©). 

- The 8.3 /im extinction as derived from the model agrees 
with the extinction derived from MSX, except for PI, 
whose lack of 8.3 fim extinction in MSX is not in accor- 
dance with the model results. The cores' clumpy nature af- 
fects their appearance in MSX but the discrepancy is not 
fully explained. Thus, mid-IR extinction alone cannot iden- 
tify all cores, especially if star formation is occurring. 
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